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Abstract 

[ Understanding the dependence of entanglement entropy on the renormalized mass in 

quantum field theories can provide insight into phenomena such as quantum phase transi- 
tions, since the mass varies in a singular way near the transition. Here we perturbatively 
calculate the entanglement entropy in interacting scalar field theory, focussing on the de- 
pendence on the field's mass. We study A0 4 and c/</> 3 theories in their ground state. By 
tracing over a half space, using the replica trick and position space Green's functions on 
the cone, we show that space-time volume divergences cancel and renormalization can be 
consistently performed in this conical geometry. We establish finite contributions to the en- 
tanglement entropy up to two-loop order, involving a finite area law. The resulting entropy 
is simple and intuitive: the free theory result in d — 3 (that we included in an earlier publi- 
cation) AS ~ 4m 2 ln(m 2 ) is altered, to leading order, by replacing the bare mass m by the 
renormalized mass m r evaluated at the renormalization scale of zero momentum. 
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1 Introduction 



Entanglement entropy is an important property of quantum mechanical states. By tracing out 
a region of a system, the resulting von Neumann or entanglement entropy S = — Trfplnp] is a 
measure of one's inability to describe the full system if one only has access to a subsystem. There 
are other ways to characterize entanglement, but the entropy is particularly fruitful. This fun- 
damental object appears in various contexts, including quantum field theory, condensed matter 
physics, black hole thermodynamics, holography, and other regimes in which quantum mechan- 
ics plays a central role. For instance, at so-called quantum phase transitions, the entanglement 
entropy can act as a diagnostic for the nature of the phase transition. Entanglement entropy 
has been the subject of various recent investigations (including Refs. [H [U El [U O [U [TJ [H]). 

A central feature of the entanglement entropy of fields in spatial dimensions d > 1 is that 
the entropy is strictly cutoff dependent (unlike the d = 1 case which involves only a logarithm 
[HI [TO]). In d = 3 it was shown in Ref. (and earlier work in |12j ) that for a scalar field in its 
ground state the entanglement entropy between the interior and exterior of a sphere is 

S ~ + subleading (1) 

where A is the surface area of the sphere and e is some microscopic cutoff on the field theory. 
(This area law can be altered to include logarithmic corrections for more complex field theories, 
such as those involving fermions |13j.) It is important to note that the entanglement entropy is 
necessarily cutoff dependent and sensitive to the detailed microphysics. For instance, in a lattice 
field theory it would depend on the lattice spacing and arrangement. Instead if one wishes to use 
the entanglement entropy to characterize UV cutoff independent properties of a system, such as 
the IR behavior at or near a phase transition, one needs to establish such sub-leading corrections 
to the entropy. In fact such contributions do exist, as has been examined in the literature (e.g., 
see [HI dSKUl IT]). 

Recently in Ref. [18] we studied free scalar field theory in D = d+1 space-time dimensions. 
The entanglement entropy was computed for a collection of geometries, including a waveguide 
with various boundary conditions. In addition to the usual cutoff dependent area piece ~ 
Aj_/e d ~ l , where A± is the d — 1 dimensional area of the transverse space and e is a microscopic 
cutoff, there is also a finite and cutoff independent area law 

A5< = / 7ci A±_ m d_1 ln(m 2 ), for d odd 
1 "1dA_\_m d ~ l , for d even 

where = (— l)( d+1 )/ 2 [12 (Air)^ 1 ^ 2 ((d — 1)/2)!] _1 for an odd number of spatial dimensions d 
and 7d = (— l) rf / 2 [12 (27r)^ _2 ^ 2 (d— l)!!] -1 for an even number of spatial dimensions d. This was 
computed in Ref. [18] using heat kernel methods, which were earlier used in Ref. [TO]. There 
was also found further sub-leading contributions, including a perimeter term and a curvature, 
or topological, term. These contributions are interesting in the sense that not only does this 
constitute a UV cutoff independent contribution that scales as the area, it also varies with 
the field's mass. In the vicinity of a phase transition, the mass varies in a singular fashion, 
and therefore so does this contribution to the entropy (by "singular" we mean that derivatives 
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diverge near the transition, though the mass itself approaches zero). So in principle, such a 
contribution can discriminate between phases. 

It is important to note that our previous results were computed only in the very special case 
of free field theory. It is natural then to ask to what extent do such results persist in interacting 
field theories. So here we would like to extend these results to the interacting case. In particular, 
we will focus on scalar field theories of the form 

C = \{dtf - imV - ^ - A*4 ( 3 ) 

i.e., renormalizable theories in d < 3, and will study them perturbatively for small g, A. We 
are especially interested in d < 3 here, so we can track all the various divergences that emerge 
in the calculation. For d > 3, the theory is non-renormalizable, and one would need to track 
new divergences. In the latter case, our results would still properly account for the leading 
divergences, but there may be sub-leading ones. 

For definiteness, consider the d = 3 interacting case. Here A is dimensionless and g has units 
of mass (H = c = 1). This implies that it would be dimensionally correct to have area law 
corrections to eq. ([!]) of the form 

S\ ~ A-~, or S\ ~ g 2 Ah\e, etc (4) 

(possibly multiplied by further dimensionless logarithmic factors), where the "1" subscript in- 
dicates that these are leading order in the couplings. If such contributions did exist, it would 
compromise the importance of the previously stated results of eq. ([2]). Recall that the signifi- 
cance of the free theory result in ([2]) is that they are manifestly cutoff independent. This means 
that as the IR parameters of the field theory vary, say in the vicinity of a phase transition, 
then the cutoff dependent piece ~ A/e 2 would drop out, leaving the purely mass dependent and 
cutoff independent contributions. On the other hand, if in the interacting theory, contributions 
such as that suggested on dimensional grounds in eq. ((4|) existed, then as one varies the IR 
parameters, such as A or g, the entropy would vary in an explicitly cutoff dependent way, due to 
the e dependence. 

In this work we show, to two-loop order, that the contributions of the form @ in fact do 
exist when written in terms of the bare parameters, in addition to the contributions in ([HE]). 
However, the main result of this work is that we show, once renormalization has been consistently 
performed, the contributions in (|J|) are absorbed into flat space renormalization of the field's 
mass. The primary change from the free theory result is that the entropy in ([2]), written in 
terms of the bare mass m, is replaced by the same result in terms of the renormalized mass m r 
evaluated at the renormalization scale of zero momentum. So as one varies the IR parameter 
m r , there is a calculable piece that changes, while the UV cutoff dependent pieces do not. This 
leaves a cutoff independent area law contribution to the entropy, even in the interacting case. 

The organization of this paper is as follows: In Section [2] we introduce the density matrix 
and partition function on the cone, as is demanded by the so-called replica trick. In Section [3] 
we study the Green's function on the cone. In Section [J] we recapitulate the free theory result. 
In Section [5] we study A</> 4 theory to leading non-trivial order in A. In Section [6] we connect the 



3 



results to mass renormalization in flat space. In Section [7] we extend these results to gcft 3 theory. 
Finally, in Section [8] we discuss our results. 



2 Density Matrix and Partition Function 

The fundamental starting point for the study of entanglement entropy is the full density matrix 
of a system p. When we only have access to a subsystem A, we focus on the reduced density 
matrix p = Tta[/5] after tracing out the complimentary region A. In this paper we shall focus 
on the basic problem where A and A represent half-spaces. Their dividing boundary is a flat 
space of dimension d± = d — 1 (counting spatial dimensions) . Although other geometries are 
of interest, this is normally the basic starting point since any smooth geometry is locally flat. 
Various physical quantities can be computed from the final trace over region A of some power 
of p. For example, the Renyi and Tsallis entropies are defined as 

In Tr \o n ] Tr \o n ] - 1 

1 — n 1 — n 

(where "Tr" means "Tr A " here). These quantities are normally computed for integer values 
n = 2,3, . . .. The trace of the nth power can be obtained by introducing a Euclideanized time 
variable, then forming a cut over the traced out region A and performing a computation on the 
nth Riemann sheet; see Fig. [JJ (left). For the n = 1 case we need to lift n to be non-integer and 
take the n — > 1 limit. Using L'Hopital's rule this may be written as 



= -Tr[plnp] (6) 



which is the so-called replica trick. S is the entanglement entropy and will be our focus, though 
other entropies (such as Renyi and Tsallis) will be readily attainable with our methods. 

By the replica trick, the nth power of the reduced density matrix can be obtained from 
defining the field theory on a cone with deficit angle 5 = 2tt(\ — n); see Fig. Q] (right). We will 
denote all quantities on the cone with subscript n, including the partition function Z n (n — V 1 
is flat space). For the field in its ground state, the nth power of the reduced density matrix is 
given by 

In Tr [p n ] = In Z n -n In Z x (7) 



and the entanglement entropy can be expressed as 

In Z n — n In Z\ 



dn 



(8) 

n— »1 



Hence our goal will be to compute the partition function Z n , which we shall do so perturbatively 
in powers of the coupling 

In Z n = In Z nfi + In Z n>1 + ... (9) 
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Figure 1: Tracing over the left half space A requires a cut along the negative real axis (left) 
and then a computation on the nth Riemann sheet in order to compute Tr[p n ] (the vertical axis 
is Euclideanized time). For n near 1 we can perform a computation on the cone (right) with 
deficit angle 5 = 2ir(l — n), (the cone's radius is sent to infinity at the end of computation.) 



At zero temperature the partition function on the cone is given by 

Z n = f V(f>e- SE ^ (10) 



where we have Euclideanized the action. For example, the free theory (written in terms of the 
bare mass m and bare field 4>) has the formal solution 

\nZ n>0 = -ilnDet [-A + m 2 ] (11) 

The derivative of this quantity, with respect to m 2 , can be readily expressed in terms of an 
integral over the Green's function at coincidence 



^—9 ln Z n,0 
om z 



-\ f d D xG n ( X , X ) (12) 



Indeed all needed quantities can be expressed as integrals over powers of the Green's function 
on the cone G n ; a quantity that we examine in the next section. 



3 Green's Function 

Previously, in the case of the free theory, the entire analysis could be performed in terms of the 
heat kernel (or the density of states) on the cone, but the interacting case will require going to 
the Green's function as it carries more information. 



So to make progress perturbatively, we need the Green's function on the cone G n . The 
Green's function satisfies the differential equation 

(-A + m 2 )G n (x,x') = 5 D (x-x') (13) 

It is important to note that since the tip of the cone breaks translational invariance, then 
G n (x, x') is a function of both x and x' separately, even though it is only be a function of the 
distance |x — x'| in flat space. Let us write our co-ordinates as x = {r, #,xj_}, where {r, 6} 
are the polar co-ordinates on the cone and xj_ are co-ordinates on the d± = d — 1 dimensional 
transverse space. The solution for the Green's function can be expressed as the following integral 
over transverse momentum and summation (e.g., see Ref. [19j ) 



2irn J (27r) d -L f-^ J q 2 + m 2 + p 2 ± 

where J is the Bessel function of the first kind, and the coefficients are dk=o = 1> dk>i = 2. 

Although the full expression for G n for arbitrary x and x' is required to analyze the full 
theory, it is somewhat unwieldy. We shall return to this in Section [7J where it will be essential. 
However, for the purposes of the free theory and the A</> 4 theory at leading order in A, it will 
suffice to study the coincidence limit x' — > x. In this limit there are simplifications to the form 
of the Green's function. Indeed by setting x' = x we can carry out the q integral, with the result 



J V ' k=0 



(15) 



where / and K are the modified Bessel function of the first and second kind, respectively. 
Following Ref. |19| we now employ the Euler-Maclaurin formula 

E«) = 2/ dk F(k)--F'(0)-2Y,j^F {2 '- 1) (0) (16) 
k=0 J ° j>l ^ 

(where B<ij are the Bernoulli numbers) to formally re-express the sum over k in terms of an 
integral over k plus a sum over a dummy index j. This leads to 

G " (X ' X) = 2^/^pT 2 J dk J */» (V m2 + p ± r ) K V» (V m% + p ^ r ) 

+ {3 > 1) (17) 



Here the j > 1 contributions to the summation do not lead to divergences at the order we are 
working, and do not concern us here. In the first term we can rescale the dummy variable of 
integration to k = k/n in order to extract the n-dependence. This allows us to re-express the 
integral in terms of G\, namely 

G n (x,x) = C? 1 (0) + / n (r) (18) 
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where we have introduced the function 



M s £^J^<fi*^') + <*>v (19) 

The first term on the right hand side in eq. (|18p is the usual flat space Green's function, which is 
formally divergent at coincident points. Since the flat space Green's function carries translational 
invar iance, it is only a function of |x — x'|, giving 

Gi(x,x) = G 1 (|x-x|) = Gi(0). (20) 

The second term in eq. f)18|) is finite and depends on radius (distance from the tip of cone) . Since 
Kq decays exponentially fast at large r, the Green's function on the cone G n rapidly approaches 
the Green's function on flat space G\ at large r, as expected. 



4 Free Theory Result 

For the free theory (or zeroth order in the couplings in the interacting theory) we can use the 
expression in eq. (|12p to determine the derivative of the density matrix as follows 

O r\ r\ 

— ^lnTr^] = lnZ n . -n— ^lnZx.o (21) 
om z om^ om* 

= -~ j d D xG n (x,x) -n J d D xG 1 (0) 



22) 



where the first integral is over the cone (as implied by the n subscript on the integral sign) 
and the second integral is over flat space; so the angular integral gives a factor of 2nn in the 
first integral and a factor of 2tt in the second integral. Substituting in (|18p . we see that the 
G\ term cancels. This is very important because such term would otherwise scale with the 
space-time volume multiplied by a UV divergent factor. Such space-time volume divergences 
are incompatible with the area law scaling and should always cancel out; a point we will return 
to later. 

After the cancellation, we are left with 

2 [■ f-oo (• id 



The integral over r can be easily performed using J °° dyy K^{y) = 1/2, and the integral over 
the transverse space is J d d± x±_ = Aj_, giving 

9 1 rp r ni 1 - ^ A f d<i± P± 1 fOA\ 

lnTr[ Po ] = -— — A ± / ... — ■ , (24) 



dm 2 11 0i 24n J (2vr) d -L m 2 + p 2 

Then by integrating this with respect to m 2 , and performing the n derivative, we are able to 
obtain the free contribution to the entropy. For the entanglement entropy we have 

1 f d d± v i 

So = -y^ ± J ln(m 2 + p\) + const (25) 
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where "const" represents a (UV divergent) quantity, independent of mass. By performing the 
integral over the transverse momenta p± in different dimensions, and extracting only the finite 
pieces, we recover the free theory results that we reported in Ref. |18| and recapitulated in the 
introduction ([2]). 



5 Quartic Interaction 



With the free theory result in hand, we can now perturbatively compute corrections from in- 
teractions. To do so we return to the full partition function. For the quartic A</> 4 theory, the 
Euclidean partition function can be expanded to order A as 



\nZ n = In J V<pe- SE[l 



In Z nfi - — 



4! 



d D x 



x 



+ ... 



In 



^,0"^ |AG„(x,x) 2 |... 



(26) 
(27) 
(28) 



where Z n $ is the partition function of the quadratic theory, and G n is the Green's function 
evaluated at coincident points; both quantities are defined on the cone, the subscript n on the 
integral indicates that it is to be performed over the cone, and the factor of 3 comes from an 
application of Wick's theorem. 

Here we need the square of the Green's function, which we write as 



G„(x,x) 2 = G!(0) 2 + 2Gi(0)/ n (r) + f n {rf. 
So the 0(A) contribution to the density matrix is 
InTrf//];] = \nZ n ,i - n\nZ Xi i 



lnZ n ,i 
_3A 

"IT 



f d D xG n (x,x) 

J n 



f-n I d D xG 1 (0) 



Q \ roo 

- — Al _2irn J drr ^(O)/^) + f n (rf 



3A 

"IF 



n 



6n 



Gi(0) 



d dl p± 



+ 2itn 



drr f n (rY 



(29) 

(30) 
(31) 

(32) 
(33) 



(again dropping the j > 1 sub- leading terms). We note that the final term here vanishes when 
we differentiate with respect to n and then set n = 1, so the only term that contributes is the 
first term. This gives the result 



d dl ?± 



1 



(34) 



Consider the d = 3 example (d± = 2). In this case, the integral is dimensionless and does not 
lead to any power law divergences. So we have Si ~ \AG\(0) (up to logarithmic corrections). 
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In 3+1-dimensions, the Green's function has a quadratic divergence at coincidence, which we can 
regulate in position space with a microscopic cutoff e, giving Gi(0) ~ 1/e 2 . Hence, S\ ~ XA/e 2 , 
as anticipated in the introduction. 

By combining this with the free theory result in (|25p . we have the following total result for 
the entropy at 0(A) 



S 



So + S% + . . . 

' d d± p± 



-1a 

12 



(2vr) c 



ln(m 2 + p\) + 



A gi(0) 
2 to 2 + 



+ const + 



(35) 
(36) 



For a chosen number of transverse dimensions d± the integral over p± can be readily performed. 
However, an understanding of the divergences in this result, coming from the Gi(0) factor and 
the high p± regime, is quite unclear when left in this form. 



6 Mass Renormalization 

In order to understand the structure of the previous result, in particular the organization of 
divergences, it is important to re-express the result in terms of renormalized couplings, rather 
than the bare couplings that we have studied so far. 

Since the UV divergences arise from arbitrarily short distance physics, and since (away from 
the tip of the cone) the field theory on the cone recovers the field theory in flat space, then we 
should recall flat space renormalization. Recall that the renormalized mass m 2 at one-loop in 
flat space Ac/> 4 theory is given in terms of the bare mass m 2 by 

ml = - 2 + ^/(0jG' 1 (p) (37) 
= m 2 + ^G l (0) (38) 

Let us conjecture that the full result for the entropy, or at least to the order we are working, 
takes the result from the free theory So in ()25[) with the bare mass replaced by the renormalized 
mass m r , i.e., 

S = -±A ± J Hm 2 r + p\) + const (39) 

By expanding this to 0(A) we indeed recover exactly So + Si, as computed in the previous 
sections. The interplay between these two contributions has a diagrammatic representation. In 
Fig. [2] (left) is the 0(A) contribution to the vacuum entanglement entropy, while Fig. [2] (right) 
is the 0(A) renormalization of the mass. 

At higher loops the mass renormalization is more involved since the renormalized mass 
acquires dependence on the renormalization scale. However, as we are studying the ground 
state, then presumably only the p r — > limit of the renormalized mass is important. This 
suggests that the leading contribution to the entropy in the full theory is to take the free theory 
result with the bare mass replaced by the renormalized mass, using p r = as the renormalization 



9 



Figure 2: Feynman diagrams for </> 4 theory at order A: two-loop vacuum contribution to entropy 
(left) and one-loop contribution to mass renormalization (right). 



scale. To test this one could either go to 0(\ 2 ) in the loop expansion in A<^ 4 theory, or consider 
an alternate theory with momentum dependence arising at the leading loop order, such as gcj) 3 
theory, as we now explore. 



7 Cubic Interaction 

We now consider a field theory with a cubic interaction of the form 



C = \{d4>? - imV - f^ 3 (40) 
In this case the partition function can be expanded to 0(g 2 ) as 

lnZ n = lnZ nfi + / d D x [ d D x' <0(x) 3 </>(x , ) 3 > o + . . . (41) 

n J n 



2!(3!) 2 

= lnZ nfi+yj^2 J d°X j d D x'[§G n {*,*!f 

+9G n (x,x)G n (x,x , )G„(x , ,x')] + ... (42) 

The second term in the integral represents a tadpole contribution. As is well known, in order to 
consistently expand around 4> = in the quantum theory, we must introduce a linear term K(j) 
to the Lagrangian (though it would be absent classically). At one-loop the bare parameter k is 
required to take the value 

K = - 9 -G l {V) (43) 

Once this is enforced, then all tadpole contributions drop out. This can be understood as a 
cancellation between the tadpole contribution to the vacuum entanglement entropy depicted in 
Fig. [3] (left) and the generation of a one-point function depicted in Fig. [3] (right). 
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Figure 3: Tadpole diagrams for (ft 3 theory at order g 2 : two- loop vacuum contribution to entropy 
(left) and one-loop contribution to the one point function (right). 



On the other hand, the first term in the integral eq. (j42[) is very important and is connected 
to the mass renormalization. The corresponding contribution to the density matrix is 

6 5 2 



InTrK] 



2!(3!) 



! d D x f d D x'G n (x,x') 3 

J n J n 



n 



d D x 



d D x' Gi(|x-x' 



(44) 



whose diagrammatic form is given in Fig. 0] (left). 

By construction the Green's function can be formally decomposed into a flat space con- 
tribution that carries translational invariance and a correction that does not. We write this 
as 

G n (x,x') = G 1 (|x-x'|) + / n (x,x') (45) 

which is a natural generalization of the decomposition we performed earlier in eq. f)18|) . Here 
we will not need to write out the full details of f n (it can be inferred from the Euler MacLaurin 
formula), but we will specify some specific properties. In particular, away from the tip of the 
cone, f n is finite as x' — > x since the divergence is captured by the flat space divergence G\. 

We insert (|4"5|) into (JUJ and expand out the cubic piece as G n = G\ + 3Gf / n + 3Gif% + f 3 . 
The first term in this expansion ~ Gf is naturally paired with the second term in (I44p . leading 
to the following difference 

6 5 2 



2!(3!) 2 

By defining y = x 



I d D xf d D x' Gi(|x-x'|) 3 -n I d D x I d D x' Gi(|x-x'|) 3 

in Jn J J 

x ; and using translational invariance, the second term is 



nV ld n yG 1 {y) 



(46) 



(47) 
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Figure 4: Feynman diagrams for (ft 3 theory at order g 2 : two- loop vacuum contribution to entropy 
(left) and one- loop contribution to mass renormalization (right). 



where V is the (flat) space-time volume. On the other hand, the first term is different since 
it involves a double integral over the cone. However, the divergent part, which arises from 
the region x' — > x, is proportional to the volume of the cone ^ nV . This leads to a natural 
cancellation between these two volume divergences. 

Instead the leading contribution to the entropy at 0(g 2 ) is from the 3G 2 f n piece, leading to 



In Tv[p n i] = J d D xj n d D x' Gi (|x - x'|) 2 /„(x, x') 



(48) 



We can readily extract the leading contribution to this from studying the x' — > x regime. This 
gives the result 

InTrK] = igl J AC,(S) ! |A/„W (49) 



2!(3!) 2 Yin \J y y > JJ (2ir) d ^ m 2 + P 2 ± 

where we have used f n {f) = ,fri( x > x )> as given earlier in eqs. (Q21CE2])- Note that the n subscript 
on the integral over y has been dropped, since the divergence structure is away from the tip of 
the cone (except for a set of measure zero). 

Then by evaluating the derivative with respect to n at n = 1 and combining with So, we 
obtain the following result for the entanglement entropy 

S = Sb + Si + ... (51) 



m 2 + p 2 ^ 



+ const + . . . (52) 



For d = 3 the factor / d D yG 1 {y) 2 gives a logarithmic divergence, leading to a contribution 
to the entropy of the form S\ ~ g 2 Alne (multiplied by a secondary logarithmic factor), as 
anticipated in the introduction. However, as in the case of \(f> theory, this result can be recast 
in the form of mass renormalization. From Fig. [4] (right) the one-loop mass renormalization in 
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g(f) theory at renormalization scale of zero momentum is 

6m2 = Y I £^ dl{p)dl{ - p) (53) 

= 9 ljd D yG x {yf (54) 

For d < 3 we can approximate m 2 = m 2 + 5m 2 , since wave- function renormalization only induces 
finite corrections at this order (see the next section for a related comment on this in the case 
of counterterms) . By substituting this into eq. (139j) and then Taylor expanding to 0(g 2 ), we 
recover the result for the entropy as given in eq. (|52p . So we have seen a consistent canceling of 
divergences, leading to eq. (|39l) . This gives the leading order result for the entropy as simply the 
free theory result with the bare mass replaced by the renormalized mass evaluated at p r = 0. 

8 Discussion 

We have shown that the UV cutoff independent area law that exists in the free theory persists 
in the interacting theory through a renormalization of the mass. This means that as we vary the 
IR parameters, the quadratically divergent contribution Sq ~ A/e 2 (for d = 3) cancels, leaving 
a UV cutoff independent area law, given in eq. (|39|) . This was determined by using the replica 
trick, requiring analysis on a conical geometry. Renormalization on the cone is non-trivial, due 
to the curvature at the tip of the cone, however the space-time volume divergences cancel in the 
entropy (this is true for both the entanglement entropy, as well as the Renyi or Tsallis entropies, 
etc) due to an interplay between two-loop vacuum contributions to the entropy and one-loop 
flat space mass renormalization. (It is common for a vacuum diagram at loop order L to be 
related to certain Appoint functions at loop order L — 1; for another scalar field example of this 
see [20]). 

The renormalization is readily seen in terms of the bare fields/couplings, as we have demon- 
strated here. However, we have also performed the calculation and obtained the same results 
with the Lagrangian re- written in terms of renormalized fields/couplings and counterterms. In 
the latter case, a new space-time volume divergence emerges in the quadratic sector due to 
wave-function renormalization, but is found to cancel against a corresponding counterterm. 

Much recent discussion of entanglement entropy has focussed on conformal field theories. 
These works are particularly interesting at a phase transition, say a quantum phase transition, 
since systems generally organize themselves into conformal field theories at such a critical point. 
Presumably, finite contributions to the entropy at a phase transition encodes some information 
about the universality class that the field theory lies in. In such cases the entropy has a divergent 
piece So ~ A/e 2 (for d = 3) and the finite pieces are 0(1) since there is no scale present, by 
definition, in order to have a finite area law. On the other hand, by introducing a mass, as 
we have done in this work, we obtain AS ~ Am 2 lnm 2 . Although this explicitly breaks the 
scale invariance, it is of interest near the phase transition. The renormalized mass m r varies in a 
singular fashion near the phase transition (singular in the sense that it's derivatives diverge), and 
hence we are extracting a singular contribution to the entropy. We have in mind the following 
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hierarchy: e 2 <C 1/m 2 , <C A, i.e., the correlation length £ ~ l/m r is much larger than the 
inter-particle spacing, as would be the case near a phase transition. But, at the same time, the 
characteristic macroscopic size of the system \[A can still be much larger than this correlation 
length. In this regime we have Am 2 . 3> 1, so the UV cutoff independent area law contribution 
to the entropy can be numerically large, and conceivably measurable [21]. 

Simple models of symmetry breaking involve a scalar field interacting under a potential 
V(<f>) = — ^a 2 + jrC 4 - Expanding around the vacuum (whether the symmetry is broken /i 2 > 
or unbroken /j? < 0) leads to a potential that fits into the form analyzed in this paper eq. ©, for 
specific values of g, m. Hence, one can apply the cutoff independent entropy law to symmetry 
breaking models, including application to Ginzburg-Landau type treatments of phase transitions, 
which would be of interest. 

It is important to extend these results, which may be possible in several directions. Firstly, 
here we have focussed on a half space geometry, though other geometries are of interest. For 
instance, one may consider a waveguide, as we studied earlier in Ref. [18], and perform some 
expansion in powers of area, perimeter, curvature, etc (perhaps similar to [221 I23j). or one may 
consider finite intervals or spheres, etc. Secondly, another important extension is to consider 
other interactions and other fields, fermions or gauge bosons (e.g., [231 [251 [26] ) . Finally, it is 
of interest to extend the results to higher loop order in the perturbative expansion, as well as 
tracking other finite sub- leading corrections. 
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